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The finite-temperature phase diagram of the attractive Hubbard model is studied by means of 
the Dynamical Mean Field Theory. We first consider the normal phase of the model by explicitly 
frustrating the superconducting ordering. In this case we obtain a first-order pairing transition 
between a metallic phase and a paired phase formed by strongly coupled incoherent pairs. The 
transition line ends in a finite temperature critical point, but a crossover between two qualitatively 
different solutions still occurs at higher temperature. Comparing the superconducting and the 
normal phase solutions, we find that the superconducting instability always occurs before the pairing 
transition in the normal phase takes place, i.e., Tc > Tpairing- Nevertheless, the high-temperature 
phase diagram at T > Tc is still characterized by a crossover from a metallic phase to a preformed 
pair phase. We characterize this crossover by computing different observables that can be used to 
identify the pseudogap region, like the spin susceptibility, the specific heat and the single-particle 
spectral function. 

PACS numbers: 71.10.Fd, 71.10.-w, 74.25.-q 



I. INTRODUCTION 

The attractive Hubbard model represents an unvalu- 
able tool to understand properties of pairing and super- 
conductivity in systems with attractive interactions. The 
simplifications introduced in this model allow a compre- 
hensive study of the evolution from the weak-couphng 
regime, where superconductivity is due to BCS pairing 
in a Fermi liquid phase, and a strong couphng regime, in 
which the system is better described in terms of bosonic 
pairs, whose condensation gives rise to superconductiv- 
ity (Bose Einstein (BE) superconductivity )i. It has been 
convincingly shown that such an evolution is a smooth 
crossover and the highest critical temperature is achieved 
in the intermediate regime where none of the limiting 
approaches is rigorously validi^. A realization of such 
a crossover scenario has been recently obtained through 
the development of experiments on the condensation of 
ultracold trapped fermionic atoms'^ . In these systems the 
strength of the attraction can be tuned by means of a tun- 
able Fano-Feschbach resonance, and the whole crossover 
can be described^. 

In the context of high-temperature superconductivity, 
the intermediate-strong coupling regime in which inco- 
herent pairs are formed well above the critical tempera- 
ture has been invoked as an interpretation of the pseu- 
dogap phased. Moreover, since the early days of the dis- 
covery of these materials, the evolution with the doping 
level of both the normal- and the superconducting-phase 
properties induced some authors^'^ to recognize the fin- 
gerprints of a crossover between a relatively standard 
BCS-like superconductivity in the overdoped materials 
and a strong-coupling superconductivity associated to 



Bose-Einstein condensation (BE) in the underdoped ma- 
terials. Indeed at optimal doping the zero-temperature 
coherence length is estimated to be around 10—20 AL^, 
i.e., much smaller than for conventional superconductors 
but still large enough to exclude the formation of local 

pairs 

It is understood that the attractive Hubbard model has 
not to be taken as a microscopic model for the cuprates, 
since a realistic description of the copper-oxygen planes of 
these materials unavoidably requires a proper treatment 
of strong Coulomb repulsion. This simplified model rep- 
resents instead an ideal framework where the evolution 
from weak to strong coupling can be studied by simply 
tuning the strength of the attraction. The main aim of 
the present work is to identify if, and to which extent, at 
least some aspects of the phenomenology of the cuprates 
can be interpreted simply in terms of a crossover from 
weak to strong coupling. 

The main simplifications introduced by the attractive 
Hubbard model can be summarized as (i) Neglect of 
repulsion. Even if some attraction has to develop at 
low energy, the large short-range Coulomb repulsion im- 
plies that the interaction must become repulsive at high- 
energy in real systems. In some sense, an attractive Hub- 
bard model picture can at most be applied to the low- 
energy quasiparticles. (ii) The model naturally presents 
s-wave superconductivity, as opposed to the d-wave sym- 
metry observed in the cuprates (Hi) Neglect of retarda- 
tion effects. The Hubbard model describes instantaneous 
interactions, while every physical pairing is expected to 
present a typical energy scale. 



2 



The model is written as 

<ij>(7 i ^ ' ^ ' 

-li'Y^{n,^ +n,{) (1) 

i 

where {cia) creates (destroys) an electron with spin 
a on the site i and Uia = cl^Cia is the number operator; 
t is the hopping amplitude and U is the Hubbard on-site 
attraction (we take U > 0, with an explicit minus sign 
in the hamiltonian) . Notice that, with this notations, 
the Hamiltonian is explicitly particle-hole symmetric for 
^ = 0, which therefore corresponds to n = 1 (half-fiUing) . 

Despite its formal simplicity, this model can be solved 
exactly only in d = 1 , while in larger dimensionality ana- 
lytical calculations are typically limited to weak {U <C t) 
or strong {U ^ t) coupling, where the BCS and the 
BE approaches are reliable approximations. It is anyway 
known that for d > 1, the ground state of Q is super- 
conducting for all values of U and all densities n, with 
the only exception of the one-dimensional half-filled case. 
At half-filling the model has an extra-symmetry and the 
superconducting and the charge-density-wave order pa- 
rameters become degenerate. 

A reliable description of the evolution of the physics as 
a function of U requires to treat the two limiting regimes 
on equal footing overcoming the drawbacks of perturba- 
tive expansions. Quantum Monte Carlo (QMC) simula- 
tions represent a valuable tool in this regard, and they 
have been applied to the twoSiiLiSiiiii and threelS di- 
mensional attractive Hubbard model. Even if the sign 
problem does not affect these simulations, finite size ef- 
fects and memory requirements still partially limit the 
potentiality of this approach. 

A different non perturbative approach is the Dynam- 
ical Mean-Field Theory (DMFT), that neglects the spa- 
tial correlations beyond the mean field level in order to 
fully retain the local quantum dynamics, and becomes 
exact in the limit of infinite dimensiongi^. Due to the 
local nature of the interaction in the attractive Hubbard 
model, we expect that the physics of local pairing is well 
described in DMFT. Moreover, this approach is not bi- 
ased toward metallic or insulating states, and it is there- 
fore particularly useful to analyze the BE-BCS crossover. 
On the other hand, the simplifications introduced by the 
DMFT are rigorously valid only in the infinite dimen- 
sionality limit, and even if the DMFT has obtained many 
successes for three dimensional systems, its relevance to 
lower dimensionality like d = 2 is much less established, 
and represents a fourth limitation of our study in light 
of a comparison with the physics of the cuprates. In 
particular, the role of dimensionality in determining the 
pseudogap properties of the attractive Hubbard model 
has been discussed in Refs^. 

The study of the attractive Hubbard model can greatly 
benefit of a mapping onto a repulsive model in a magnetic 
field. The mapping is realized in a bipartite lattice^® by 



a 'staggered' particle-hole transformation on the down 
spins Cii (— l)*c||. The attractive model with a finite 
density n transforms into a half-filled repulsive model 
with a finite magnetization m = n — 1. The chemical po- 
tential is transformed, accordingly, into a magnetic field 
h = ji. In the n = 1 case (half-filling) the two models 
are therefore completely equivalent. We notice that the 
above mapping does not only hold for the normal phases, 
but extends to the broken symmetry solutions. The three 
components of the antiferromagnetic order parameter of 
the repulsive Hubbard model are in fact mapped onto a 
staggered charge-density-wave parameter (z component 
of the spin) and an s-wave superconducting order param- 
eter {x—y components). The above mapping is extremely 
useful, since it allows to exploit all the known results for 
the repulsive model and for the Mott- Hubbard transition 
to improve our understanding of the attractive model. 

In recent works the DMFT has been used to study the 
normal phases of the attractive Hubbard model. In par- 
ticular, a phase transition has been found both at finiteiS 
and at zero temperature^*' between a metallic solution 
and an pairing phase of pairs. The insulating pairs phase 
is nothing but a realization of a superconductor without 
phase coherence, i.e., a collection of independent pairs. 
As it has been discussed in Ref.^^'^", this phase is the 
'negative-C/' counterpart of the paramagnetic Mott insu- 
lator found for the repulsive Hubbard model. We notice 
that the insulating character of the pairing phase is a 
limitation of the DMFT approach, in which the residual 
kinetic energy of the preformed pairs is not described. 
The pairing transition has been first identified in RefiiS 
by means of a finite temperature QMC solution of the 
DMFT. The T = study of Ref.^o has clarified that the 
pairing transition is always of first order except for the 
half-filled case, and that it takes place with a finite value 
of the quasiparticle weight Z = (1 — d'S{uj) / duj)~^ , asso- 
ciated to a finite spectral weight at the Fermi level. In 
the latter paper, it has also been shown that the pairing 
transition gives rise to phase separation. 

For what concerns the onset of superconductivity, a 
DMFT calculation of the critical temperature Tc has 
been performed for the case of n = 0.5 in the same 
Refii^. The Tc curve, extracted from the divergence of 
the pair-correlation function in the normal phase, dis- 
plays a clear maximum at intermediate coupling and re- 
produces correctly both the BCS and the BE predictions 
in the asymptotic limits, remaining finite for all U =/= 0. 

In this work we complement the analysis of Ref.^-, 
by extending our phase diagram to finite temperature, 
still using Exact Diagonalization (ED) to solve the im- 
purity model associated with the DMFT of the Hubbard 
model^i. We also compare the normal state solutions 
with the superconducting solutions which are stable at 
low temperatures. The use of ED allows us to reach ar- 
bitrarily small temperatures which are hardly accessible 
by means of QMC. Quite naturally, the extension of ED 
to finite temperature requires a more severe truncation of 
the Hilbert space. We have checked that all the thermo- 
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dynamical quantities we show are only weakly dependent 
on the truncation. The plan of the paper is the following: 
in Sec. II we briefly introduce the DMFT method and its 
generalization to the superconducting phase; In Sec. Ill 
we discuss the finite temperature phase diagram in the 
normal phase characterizing the low-tcmperature pairing 
transition; In Sec. IV we analyze the superconducting so- 
lutions; In Sec. V we compare different estimators of the 
pseudogap temperature in the high-temperature normal 
phase. Sec. VI contains our concluding remarks. 

II. METHOD 

The DMFT extends the concept of classical mean-field 
theories to quantum problems, by describing a lattice 
model in terms of an effective dynamical local theory. 
The latter can be represented through an impurity model 
subject to a self-consistency condition, which contains 
all the information about the original lattice structure 
through the non-interacting density of states (DOS)ii&. 
Starting from the Hubbard model we obtain an at- 
tractive Anderson impurity model 

Ti-AM = - ^ Vfe4^^co,cr + H.c. + ^ efc4_^Cfc,o. 

k,a k.(7 

The self-consistency is expressed by requiring the identity 
between the local self-energy of the lattice model and the 
impurity self-energy 

S(ic^„) = g°(zc^„)-i-G(ic^„)-\ (3) 

where G(iw„) is the local Green's function of and 
5°(ia;„)~^ is the dynamical Weiss field, related to the 
parameters in |2Jl by 

g°(ic^„)-^ ^zc^„ + M- V . . (4) 

By expressing the local component of the Green's func- 
tion in terms of the lattice Green's function, namely 
G{r — 0, iujn) — J2k '^i^' i^n), Eq. Q implies 

g\^u,r' = ( I de . ^^"^ \ ' + E(*c.„), 

(5) 

where D{e) is the non interacting density of states of 
the original lattice. We consider the infinite-coordination 
Bethe lattice, with semicircular DOS of half-bandwidth 
D (i.e., D{e) = {2/7rD^)VD^ - e^), for which Eq. © is 
greatly simplified and becomes 

g°{iuin) ^=icJ„ + M -^Giiun)- (6) 

In this work we also consider solutions with explicit s- 
wave superconducting order, by allowing for local anoma- 
lous Green's functions F{t) = — (TtCo|(t)cox). The 



whole DMFT formalism can then be recast in Nambu- 
Gorkov spinorial representationi&, and Eqs. Q and Q 
must be read as matrix identities in the Nambu space. 
As far as the impurity model is concerned, we need to de- 
scribe an Anderson impurity model with a superconduct- 
ing bath or, equivalently, with an anomalous hybridiza- 
tion in which Cooper pairs are created and destroyed in 
the electronic bath, i.e., a term V^f (cfe|Cfcj -I- H.c.) is 
added to (0). 

The heaviest step of the DMFT approach is to com- 
pute G{iLUn) for the Anderson model This solution 
requires either a numerical approach or some approxi- 
mation. Here we use Exact Diagonalization. Namely, we 
discretize the Anderson model, by truncating the sums 
over k in Eqs. lU and Q to a finite number of levels 
Ns . It has been shown that extremely small values of Ns 
provide really good results for thermodynamic properties 
and reliable results for spectral functions. In this work 
we use the ED approach at finite temperature, where it is 
not possible to use the Lanczos algorithm, which allows 
to find the groundstate of extremely large matrices. To 
obtain the full spectrum of the Hamiltonian, needed to 
compute the finite-temperature properties, we are forced 
to a rather small value of Ns, up to 6. All the results 
presented here are for Ns = 6, and we always checked 
that changing Ng from 5 to 6 does not affect the relevant 
observables we discuss in the present work, except for the 
real-frequency spectral properties. 



III. THE PAIRING TRANSITION 

In this section we limit our analysis to normal phase 
paramagnetic solutions in which no superconducting or- 
dering is allowed. Even if the s-wave superconducting 
solution is expected to be the stable one at low tempera- 
tures, our normal state solutions are representative of the 
normal phase above the critical temperature. The region 
in which the normal state is stable may of course be en- 
larged by frustrating superconductivity through, e.g., a 
magnetic field. Moreover, the nature of the normal phase 
gives important indications on the nature of the pairing 
in the different regions of the phase diagram. As men- 
tioned above, it has been shown that the normal phase 
of the attractive Hubbard model is characterized by a 
"pairing" transition between a Fermi-liquid phase and a 
phase in which the electrons are paired, but without any 
phase coherence among the pairs. 

The pairing transition has been first discussed at finite 
temperature in Ref.^^, and a complete characterization at 
T — has been given in RefiSS. In this paper we complete 
the finite temperature study of the transition and connect 
it to the zero-temperature phase diagram, finally drawing 
a complete phase diagram in the attraction-temperature 
plane for a density n = 0.75, taken as representative of 
a generic density (except for the peculiar particle-hole 
symmetric n = I case). This situation would correspond 
to a repulsive model at half-filling in an external mag- 
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netic field tuned to give a finite magnetization m = 0.25. 
The T = DMFT solution of the attractive Hubbard 
model is characterized by the existence of two distinct 
solutions, a metallic one with a finite spectral weight 
at the Fermi level and an insulating solution formed by 
pairs, with no weight at the Fermi level. The previ- 
ous study has also clarified that the quasiparticle weight 
Z = {1 — dY.{u!) / duj)^^ , which may be used as a sort 
of order parameter for the Mott transition at half-filling, 
loses this role for the doped attractive Hubbard model, 
being it finite both in the metallic and pairing phases. 
At T = 0, the metallic solution exists only for U < Uc2, 
and the insulating one for U > Uci, with Uci < Uc2- In 
other words, a coexistence region is present where both 
solutions exist, and where the actual ground state is de- 
termined minimizing the internal energy. The clear- 
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FIG. 1: Evolution of the imaginary part of the Green's func- 
tion as a function of temperature for U/D — 2.4. In each 
panel are shown the metallic and insulating ( x ) solutions, 
(the chosen value of the attraction lies in the coexistence re- 
gion). 

cut r = characterization of the two solutions based 
by the low-energy spectral weight is lost at finite tem- 
perature, where both solution have finite weight at the 
Fermi level. Nonetheless, two families of solutions can 
still be defined, each family being obtained by continu- 
ous evolution of the different T = phases. The two 
solutions are still clearly identified at relatively low tem- 
peratures, further increasing the temperature, the dif- 
ferences between the two solutions is gradually washed 
out, as shown in Fig. ^ where we plot the tempera- 
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FIG. 2: Average double occupation as a function of the at- 
traction strength for different temperatures. The first order 
transition at low temperatures becomes a continuous evolu- 
tion at high temperatures, where there is no more distinction 
between metallic and pairing solutions. 



ture evolution of the imaginary part of the Green's func- 
tion in imaginary frequency for U = 2. AD, which lies 
in the T = coexistence region. While at T = 1/75D 
and T — 1/50D the difference in the two solutions is 
still clear, at T — 1/31D, the two solutions become ba- 
sically indistinguishable. This result suggests that, as 
intuitively expected, the temperature reduces the differ- 
ence between the solutions and consequently, the size of 
the coexistence region, which is expected to close at some 
finite temperature critical point (the attractive counter- 
part of the endpoint of the line of metal-insulator tran- 
sitions in the repulsive mode]Mt). A similar information 
is carried by the analysis of the average value of double 
occupancy ~ (n^ni). This quantity naturally discrim- 
inates between an pairing phase with a large value of rid 
and a metal with a smaller value. As shown in Fig. |21 at 
low temperature we have two solutions with a different 
value of in the coexistence region, and a jump in this 
quantity at the transition. Upon increasing the temper- 
ature, the two solutions tend to join smoothly one onto 
the other, signaling again the closure of the coexistence 
region, which is substituted by a crossover region. Anal- 
ogous behavior is displayed by the quasiparticle weight 
Z. 
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Repeating the same analysis for a wide range of cou- 
pling constants and temperatures, we are able to con- 
struct a finite-temperature phase diagram for the pairing 
transition, shown in Fig. 13 For temperatures smaller 
than a critical temperature Tpairing, we compute the fi- 
nite temperature extensions of Uci and Uc2, which mark 
the boundary of the coexistence region. The two lines 
(depicted as dashed lines in Fig. 3) converge into a finite 
temperature critical point at U = Upairing — 2. 313 and 
T = Tpairing — 0.03D. Despitc the closure of the coexis- 
tence region, a qualitative difference between weak cou- 
pling and strong coupling solutions can still be identified 
for T > Tpairing , determining a crossover region in which 
the character of the solution smoothly evolves from one 
limit to the other as the attraction is tuned. At this stage, 
the crossover region is "negatively" defined as the range 
in which the Green's function does not resemble any of 
the two low temperature phases. The crossover lines are 
estimated as the points in which it becomes impossible 
to infer from the Matsubara frequency Green's function 
whether the low-energy behavior is metallic or insulating, 
ft has been shown for the repulsive Hubbard model that 
this kind of crossover is accompanied by a qualitative dif- 
ference in transport properties. In the region on the left 
of the crossover, the conduction is metallic and the re- 
sistivity increases with temperature. In the intermediate 
crossover region the system behaves like a semiconduc- 
tor with a resistivity which decreases upon heating, and 
finally in the phase on the right of the crossover region 
the system behaves like a heated insulator 

Coming from the left, the first crossover occurs when 
G{iLUn) has no longer a clear metallic behavior with a 
finite value at zero frequency, while the second crossover 
line delimits the region in which the gap of the paired 
solution is closed by thermal excitations. We will come 
back later to the crossover region and compare the above 
defined lines with physically sensible estimators of the 
pseudogap temperature, like the specific heat and the 
spin susceptibility. 

Turning to the coexistence region, we can also ask our- 
selves which is the stable phase. This requires a compari- 
son between the Gibbs free energies of the two phases. At 
half-filling, where the attractive and the repulsive model 
are equivalent, it has been shown that at T = the metal- 
lic solution is stable in the whole coexistence region^i. At 
finite temperature it has been shown numerically that the 
insulator becomes stable in a large portion of the coex- 
istence region due to its large entropj*i&. The transition 
is therefore of first order for all temperature below the 
critical temperature, except for the two second-order end- 
points at r = and T — Tpairing- For densities out of 
half-filling it has been shown in Ref. i2Q| that the transi- 
tion is of first order already at T = and it is accompa- 
nied by a small phase separation region. For n = 0.75, 
the T = first-order transition occurs quite close to Uc2- 

Analogously to the half-filling case, the finite temper- 
ature almost immediately favors the pairing phase. In- 
deed, computing the free energy following, e.g., Refn^i, 
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FIG. 3: Phase diagram in the U-T plane. At low temperature 
two critical lines Uci{T) and Uc2{T) individuate the coexis- 
tence region. The two lines converge in a finite temperature 
critical point. At higher temperatures we can still define two 
crossover lines. The superconducting critical temperature is 
also drawn as a solid line (cfr. Fig. 4). 

we find the pairing phase stable for almost every point 
in the coexistence region. We had to use an extremely 
dense mesh of points in the U direction to identify a 
small section where the metallic phase is stable at finite 
temperature. Therefore the finite temperature first-order 
transition occurs extremely close to the Ud line for finite 
temperature and rapidly moves closer to Uc2 only at re- 
ally small temperatures. 

IV. THE SUPERCONDUCTING PHASE 

The above stability analysis has been restricted to nor- 
mal phase solutions. Indeed the superconducting solu- 
tion is expected to be the stable one at T = for all 
densities and values of the interaction U. The critical 
temperature Tc is obtained directly as the highest tem- 
perature for which a non-vanishing anomalous Green's 
function F{u;) exists. 

The DMFT critical temperature Tc for n = 0.75 as 
a function of U is reported in Fig. 21 (full dots) and it 
qualitatively reproduces the limiting behavior, with an 
exponential BCS-like behavior for small U's and a 1/U 
decrease at large U according to the expression for the 
BE condensation temperature hard-core boson systemS^. 
As a result, Tc assumes its maximum value of about O.I_D 
for an intermediate coupling strength Umax — 2.1D. In- 
terestingly, the maximum Tc occurs almost exactly at the 
coupling for which the pairing transition in the normal 
phase would take place in the absence of superconductiv- 
ity. 

It might be noticed however that, while the BE re- 
sult (open triangles in Fig. basically falls on top of 
the DMFT results, the BCS formula (open circles) only 
qualitatively follows the full solution. This "asymmetry" 
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in recovering the BCS behavior arises from the partial 
screening of the bare attraction due to second order po- 
larization termsS^. Because of these corrections the at- 
traction is renormalized as Ue/f ~ — AU'^/t, so that 
~ jj{l + AU/t) = 1/U + A/t. When this correction 
is plugged in the BCS formula for Tc, it results in a correc- 
tion to the prefactor. If we simply extract the rescaling 
factor for a given small value of U {Tc/T^^'^ ~ 0.32) and 
we simply scale the whole weak-coupling curve by this 
factor, we obtain the points marked with asterisks, whose 
agreement with the DMFT results does not require fur- 
ther comments. It is interesting, instead, to compare the 
DMFT estimations for Tc with the QMC results: despite 
the presence of many factors (such as the exact shape of 
the D.O.S. of the model or the finite dimension effects) 
which are capable to introduce relevant variations in the 
values of Tc, some general similarities appear clearly. In- 
deed, while simple rescaling the data in terms of the half- 
bandwidth D, both Tc and Umax estimations with the 
twoiSii^ and three^^ dimensional QMC are lower than 
the DMFT evaluation (i.e., Tc - 0.04D, Umax ^ 0-7D 
for the d = 2 case, and Tc ~ 0.05D, Umax 1.3D in 
d = 3, even if for a lower density of n = 0.5), one can 
observe, quite surprisingly, that the ratio between Tc and 
Umax is around 0.04 0.05 in both the DMFT and the 
two QMC cases. 
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FIG. 4: Critical temperature as a function of U at n — 0.75: 
the DMFT data (black circle) are compared with the BCS 
-both the bare (empty circles) and the rescaled one (stars)- 
and the BE mean field predictions (em pty triangles) for an 
hard-core boson systems (see Refs. 

Coming back to our DMFT results, the simplest and 
most important observation is that the critical temper- 
ature is always higher than the critical temperature for 
the pairing transition. For example in our n — 0.75 case, 
T^""^ is about O.ID, against a Tpa^„ng of 0.03D. As a 
result, the whole phase transition is hidden by supercon- 
ductivity, which remains as the only real instability of the 
system (cfr. Fig. 3). Nevertheless, the crossover lines at 
higher temperature survive the onset of superconductiv- 



ity. Therefore the normal phase we reach for T > Tc is 
really different according to the regime of coupling we are 
in. For weak-coupling, the normal phase is substantially 
a regular Fermi-liquid and superconductivity occurs as 
the standard BCS instability. In the strong-coupling 
regime, the normal phase is instead a more correlated 
phase which presents a pseudogap in the spectrum. At 
intermediate coupling, where the superconducting criti- 
cal temperature reaches its maximum, the normal phase 
is in a crossover region between the two limiting behav- 
iors. 



V. THE PSEUDOGAP PHASE: SPIN 
SUSCEPTIBILITY, SPECIFIC HEAT AND 
SPECTRAL FUNCTIONS 

Even if the onset of superconductivity completely 
hides the pairing transition, the fingerprints of the low- 
temperature normal phase are still visible in the high- 
temperature phase diagram, in which a crossover from a 
metallic phase to a gapped phase is still present. It is 
tempting to associate the region in which the system be- 
haves as a collection of incoherent pairs to the pseudogap 
regime of the cuprates. It is important to underline that, 
in this framework, the definition of the pseudogap phase 
is somewhat tricky, and it implies a certain degree of ar- 
bitrariness. In this section we come back to this region 
and compute various observable whose anomalies have 
been used to identify the pseudogap phase and compare 
the related estimates of the pseudogap temperature T*. 

Our first estimate is based on the evaluation of the 
uniform spin susceptibility Xs as a function of tempera- 
ture for different attraction strengths. The opening of a 
gap in the spin excitation spectrum, not associated with 
any long-range order, represented in fact one of the first 
indications of existence of the pseudogap phase in high- 
temperature superconductors. The DMFT calculation of 
Xs can be performed by evaluating the derivative of the 
magnetization m with respect to a uniform magnetic field 
in the limit of vanishing h. In terms of the local Green 
functions 

l„Ea;„[G'T(^'^n)-G'i(ia;„)] 
- tl 2^ h ■ 

This calculations has been performed by varying the tem- 
perature in a wide range (0 < T < 2D) for four different 
values of the pairing interaction {U /D ~ 0.8, 1.8, 2.4 and 
3.6) and represent an extension of the results reported in 
Ref. . The results of our calculation are summarized 
in Fig. El 

In the weak-coupling side {U — Q.8D) we find a con- 
ventional metallic behavior of Xs, which increases mono- 
tonically with decreasing temperature. The interaction 
reduces the zero-temperature extrapolated value with re- 
spect to the non-interacting result Xs = p(0). On the op- 
posite side of the phase diagram, in the strong-coupling 
regime {U = 3.61?) the standard high-temperature be- 
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FIG. 5: Spin susceptibility in the normal phase as a func- 
tion of temperature for U/D = 0.0,0.8,1.8,2.4,3.6. These 
values of U are representative of all the interesting region of 
the phase diagram in Fig. |3] moving from the metallic to 
the paired side. The values of the superconducting critical 
temperature are marked by small black arrows. 



havior of Xs extends only down to a certain tempera- 
ture T^.j^ where a maximum of Xs is reached. When 
the temperature is further reduced Xs starts to decrease, 
exponentially approaching zero in the T — Q limit, sig- 
naling the opening of a gap in the spin-excitation spec- 
trum. A qualitatively similar behavior is found also in 
the intermediate-coupling regime, at least as long as the 
value of U stays larger than Upairing (e.g., U = 2 AD), 
or, in other words, as long as the left line defining the 
crossover region in Fig. Olis not crossed. The behavior of 
Xs becomes richer for U — 1.8-D < Upairing- At high 
temperatures Xs closely resembles the insulating case, 
displaying a clear maximum at a temperature T^.j. By 
approaching T = 0, Xs no longer vanishes, but it rises 
at small temperatures displaying a minimum for a tem- 
perature lower than T^: a metallic behavior is there- 
fore recovered, associated to the narrow resonance at the 
Fermi leveliS. Such a behavior naturally defines a dif- 
ferent temperature scale Tj^, which is associated to the 
minimum of Xs and represents the lower border of the 
pairing zone, or in a sense, of the "pseudogap" region. 
Conversely, this low-temperature behavior has not been 
observed to our knowledge in finite-dimensional QMC 
simulationSiiiiiSiiiiiii^. In practice, the system displays 
a pseudogap behavior in the region between T*^^ and T^^, 
whose boundary, labeled as T*, is represented in Fig. 8. 

We finally mention that the temperature T^j for which 
Xs is maximum scales with U . This finding is in a quali- 
tative agreement with a QMC simulationliil^, where the 
T^j{U) is taken as a definition of the temperature below 
which the pseudogap appears. From a more quantitative 
point of view, as happens for Tc and Umax, the values 
of Tlf{U) of the QMC simulations are lower than our 
DMFT results (i.e., Tl^{Umax) ~ 0.15D when d = 2 



and ~ 0.45D for d = 3, against the DMFT estimate of 
~ 0.7D) . However, also in this case the ratio between 
Tlj{Umax) and Umax has a more universal value around 
0.2^0.3. 

Another relevant quantity is the specific heat Cy = 
dE/dT = -Td^F/dT^, that we obtain by differentiat- 
ing a fit to the DMFT internal energy E{T) for the same 
attraction strengths and report in Fig. |H1 Also for this 
quantity the weak coupling case {U/D = 0.8) behaves as 
a regular metal, with a linear behavior at small temper- 
atures (Cy = with 7 cx 1/m*, m* being the effective 
mass), followed by a smooth decrease when the temper- 
ature exceeds the typical electronic energy scale. The 
same qualitative result is found for the noninteracting 
system, where the low-T slope is smaller since the inter- 
acting system has a larger effective mass. In the opposite 
strong coupling limit we observe the typical activated be- 
havior of gapped systems for small temperatures, with an 
exponential dependence of Cv{T) which extends up to a 
temperature large enough to wipe out the effect of 
the gap. It is therefore natural to associate such a tem- 
perature to the closure of the pseudogap. 

In the most interesting U — 1.8-D case, two fea- 
tures are clearly present in the Cv{T) curve. The first, 
low-temperature feature is the evolution of the small-[/ 
metallic feature, which acquires a larger slope as U/D is 
increased due to the enhancement of the effective mass, 
and shrinks as a consequence of the reduced coherence 
temperature of the metal. The second feature is instead 
the evolution of the large-t/ insulating one, and would 
show an activated behavior partially hidden by the low- 
T metallic peak. Thus, the system behaves like a metal in 
the small temperature range, while it has a pseudogap for 
intermediate temperature. We estimate the lower bound- 
ary of the pseudogap region in this intermediate coupling 
regime through the maximum of the low-temperature fea- 
ture, which is controlled by the effective coherence scale 
of the metal. The upper bound is naturally defined as 
the temperature in which the activated behavior disap- 
pears. As a result, the specific heat analysis determines 
a pseudogap region with a very similar shape than the 
one determined through the spin susceptibility, with a 
re-entrance of metallic behavior in the intermediate cou- 
pling regime at low temperatures. 

An inspection to the spectral function can strengthen 
our insight on the pseudogap phase. In principle the 
ED algorithm allows to directly compute finite frequency 
spectral functions p{lo) — —l/T:ImG{uj) , avoiding the 
problems and ambiguities intrinsic with analytic contin- 
uation techniques. Unfortunately, the discretization of 
the Hilbert space which allows for an ED solution re- 
sults in "spiky" local spectral functions formed by a col- 
lection of (5- functions. In this light, wc find it useful 
to compute p{uj) by analytically continuing Eq.©, an- 
alytically computing the local retarded Green function 
Gioc(w) = -I/tt / deD{e){LJ - e + p - S,et(cc))-i. This 
procedure provides more "realistic" descriptions of both 
the non-interacting DOS and of the strong-coupling pair- 
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FIG. 6: Specific heat as a function of temperature for U /D = 
0.8,1.8,2.4,3.6. All the Cv lines are obtained by differen- 
tiating the the internal energy Ei„t{T). The expression of 
Ei„t(T] is comDuted directly bv fitting the DMFT data. 




FIG. 7: Here are plotted the density of states p{lo) for three 
different values of the interaction: U = 0.8D, 2AD and 3.6 
(from the upper to the lower row) both at the low-temperature 
(T = O.ID, left panels) and at the high-temperature (T = ID, 
right panels). 



ing phase. However, even though the spectral functions 
are smoothed by this procedure, we can only extract in- 
formations about the gross features of the spectra as, 
e.g., the amplitude of the gap A(T). Keeping these hm- 
itations in mind, some results for p{lo) are plotted in 
Fig. [3 for the weak- {U = 0.81?), the intermediate- 
[U — 2AD) and strong-coupling {U = i.QD) case a low 
and high-temperature set of data are shown. Apart from 
the obvious appearing and enlarging of a gap in p{uj) with 
increasing J7, which is evident in the low-temperature 
data, it should be noticed that both in the intermedi- 
ate and the strong-coupling regime there is apparently 
no tendency to a 'closure' of such a gap when the tem- 
perature is raised. Indeed, as it is shown in the sec- 



ond and the third row of Fig. [7| the gap starts to fill 
at some temperature (T ~ 0.451? for U = 2AD and 
T - 1.5 ^ 2.0Z? for U = 3.6D), but for these values of U 
much of the spectral weight remains in the high-energy 
Hubbard bands, and the gapped structure does not com- 
pletely vanish up to the highest temperature reached in 
our calculation [T ~ 2D). On the other hand, QMC 
results in d = 2^ obtained through maximum entropy 
show a closure of the gap in p{ui) at a temperature lower 
that our threshold. Further investigation is needed to 
understand whether the discrepancy is due to a different 
behavior between d — 2 and the infinite dimensional- 
ity limit, or it is determined by the technical difficulties 
involved in the calculation of real frequency spectra in 
both approaches. The persistence of the gap structure 
at high temperature that we find in DMFT is also ob- 
tained within a perturbative analysis of superconducting 
fluctuations at strong coupling in d = 3^. 

In Fig. |Slwe compare our estimates of the pseudogap 
temperature obtained through different physical quan- 
tities. We draw the borders of the pseudogap region 
as determined from the spin susceptibility {T*{U)) and 
the specific-heat behavior {T^{U) and the value of the 
superconducting gap Aq at zero temperature: The up- 
per borders of the spin and the specific-heat "pseudo- 
gap" region scale roughly with C/, as Aq does, so that 
both Tl{U) and T*{U) are proportional to Aq, as the 
experimentally determined pseudogap. At low tempera- 
ture, the pseudogap region boundary as extracted from 
thermodynamic response functions displays a clear re- 
entrance, which can be associated with the onset of the 
low-temperature quasi-particle peak. We also notice that 
the low-temperature curve qualitatively follows the be- 
havior of the Uc2{T) line. As mentioned above, the slope 
of Uc2 (T) is easily interpreted in terms of entropy balance 
between the two phases, which favors the preformed pairs 
phase. 

Our phase diagram also represents a warning regarding 
attempts to extrapolate the low-temperature behavior 
from the high-temperature data in order to compare with 
finite-dimensional QMC calculations. If one, as, e.g., 
in Ref. 15j extrapolated the high-temperature behavior 
down to T = in order to estimate the metal-insulator 
point, would have obtained an estimate of U* signifi- 
cantly lower than the real Uc2- This finding emphasize 
how the high-temperature properties of the attractive 
Hubbard model are only weakly dependent on dimension- 
ality, as indicated by the similarity between DMFT and 
finite-dimension QMC, while the low-temperature behav- 
ior may well be dependent on the dimensionality, as well 
as on the details of the bandstructure of the underlying 
lattice. 



VI. CONCLUSIONS 

In this paper we have investigated the finite tempera- 
ture aspects of pairing and superconductivity in the at- 
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FIG. 8: Different estimates of the pseudogap temperature. 
The filled circles indicates Ts{U), i.e., the temperatures of 
both the maxima and the minima of Xs(T'), while the empty 
circle marked ThiU) that is the temperature associated with 
the maxima and the minima of CviT). The regions on the 
left of these two lines can be interpreted as the zone of "pseu- 
dogap" behavior for the spin and the specific heat respec- 
tively. These lines are then compared with the behavior of 
the anomalous part of the self-energy at zero temperature 
(Ao, empty triangles). 



close to the zero-temperature critical point of the normal 
phase. 

The presence of the pairing transition affects the nor- 
mal phase above Tc also when superconductivity estab- 
lishes. In particular, one could be tempted to identify 
the phase of preformed pairs obtained at strong cou- 
pling with the pseudogap behavior observed in cuprates. 
In order to test the adequacy of such an identification, 
we computed different observables, whose anomalies can 
identify the appearance of the pseudogap, like the spin 
susceptibility, the specific heat and the single particle 
spectral functions. In the intermediate region of cou- 
pling, where the pairing transition occurs and the su- 
perconducting critical temperature reaches its maximum, 
the pseudogap region presents a re-entrance at low tem- 
peratures associated with a small coherent peak in the 
spectral function. At temperatures smaller than this co- 
herence temperature the system behaves like a normal 
metal with renormalized effective mass. On the other 
hand, the high-temperature boundary of the pseudogap 
region scales with U regardless the criterion we use to 
estimate it. The estimate of the pseudogap temperature 
from specific heat and spin susceptibility both scale with 
the zero-temperature gap, as in the cuprates. 



tractive Hubbard model by means of DMFT, considering 
both normal and superconducting solutions. 

In the normal phase we have identified two families of 
solutions, a Fermi-liquid metallic phase and a preformed- 
pair phase with insulating character. The latter phase is 
formed by local pairs without phase coherence. A fi- 
nite region of the coupling-temperature phase diagram is 
characterized by the simultaneous presence of both solu- 
tions. In the low temperature regime a first-order tran- 
sition occurs within this region when the free energies of 
the two solutions cross, and the region closes at a cer- 
tain temperature {Tpairing — 0.031?) in a critical point. 
Interestingly, some trace of the two solutions survives 
even for temperature larger than the critical tempera- 
ture, and two crossover lines can be defined separating a 
normal metal, a sort of semiconductor in which the gap 
is closed by temperature, and the preformed-pair phase 
with a well defined gap. 

When superconductivity is allowed, the superconduct- 
ing solution is stable for all values of the attraction and 
the critical temperature is always larger than the pair- 
ing transition temperature in the normal phase. In the 
superconducting state, we find an evolution from a weak- 
coupling BCS-like behavior, with exponentially small Tc 
from a normal metal to the superconductor, and a strong- 
coupling regime in which superconductivity is associated 
to the onset of the phase-coherence among the preformed 
pairs that occurs at Tc oc /U . The highest Tc is ob- 
tained in the intermediate region between this two lim- 
iting cases, namely for U ~ 2. 11?, which is extremely 



The most striking difference between our pseudogap 
phase-diagram and the experiments in the cuprates is 
that the pseudogap phase in the attractive Hubbard 
model is much larger than the experimental one, as it 
is measured by the large value of T*^lTc ~ 5 at the opti- 
mal value of the attraction. The experimental T* around 
optimal doping is instead very close to Tc, and, according 
to some authors, the pseudogap line tends to zero at op- 
timal doping. Moreover, the pseudogap temperature ob- 
served in the cuprates is definitely much smaller than the 
one found within our DMFT of the attractive Hubbard 
model. This inadequacy of the attractive Hubbard model 
in describing some features of the pseudogap phase de- 
scend from the above mentioned strong simplifications of 
the model (neglect of retardation effects. Coulomb repul- 
sion and d-wave symmetry of the gap) and of our DMFT 
treatment which is exact only in the infinite dimension- 
ality limit. One could be tempted in maintaining an at- 
tractive Hubbard model description for the quasiparticles 
alone, but it is important to point out that this interpre- 
tation can not be pushed too far. As an example, it is 
clear that such a description would fail for temperatures 
larger than the quasiparticle renormalized bandwidth. 

A better description of the pseudogap phase would re- 
quire models in which both an attraction and a repulsion 
are present. This is for instance the case of the mod- 
els introduced in Refs. I23i2^ where the superconducting 
phenomenon only involves heavy quasiparticles which ex- 
perience an unscreened attraction and a richer behavior 
of the pseudogap (which in this case closes around opti- 
mal doping) is found. 
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